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ABSTRACT 


An analog computer study is made of the low-frequency flow dynamics observed with 
liquid hydrogen flowing in two different nuclear -rocket cold-flow (unfueled) engine sys- 
tems. The prediction of the low-frequency response depends largely upon the generation 
of the proper mass-flow storage in the system. This is achieved by relating, in the 
model lumping, the average values of fluid enthalpy and a certain thermodynamic fluid 
property. A two -phase state equation is derived for thermodynamic equilibrium between 
the gas and liquid phases. Parametric studies are made which demonstrate the effect 
of the pump characteristics on the flow damping. 
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AN ANALOG COMPUTER STUDY OF THE LOW-FREQUENCY FLOW DYNAMICS 
OF TWO NUCLEAR-ROCKET COLD-FLOW ENGINE SYSTEMS 
by Harry W. Fox. Jr., and Ronald J. Blaha 
Lewis Research Center 

SUMMARY 


A study is made of the mechanisms which govern the low-frequency flow dynamics 
observed with liquid hydrogen flowing in two nuclear-rocket cold-flow (unfueled) engine 
systems. An analog computer model is used which is characterized by three fluid-lumps: 
a liquid lump, a two-phase lump, and a gas lump. The positions of the fluid-lump 
boundaries are calculated in high-speed repetitive operation and are allowed to shift in 
time. 

A two -phase equation of state is derived for the case when thermodynamic equilib- 
rium may be assumed between the gas and liquid phases. The form of the equation is 

P 2p = p 2p[“< P 2p> * h 2p) 


where Pg p is the two-phase pressure, is the two -phase density, 1!^ is the two- 
phase enthalpy, and a and /3 are coefficients which are functions of pressure. 

An important result of the study is that prediction of the low-frequency response 
(<5 Hz) depends largely on the generation of the proper mass-flow storage in the system. 
It is found that most of the mass-flow storage occurs in two-phase regions of very low 
mass quality. A good representation of this mass-flow accumulation, in the model 
lumping, is achieved by relating the average fluid enthalpy to the average value of the 
thermodynamic property [dp/d P]^ which is the partial derivative of density with respect 
to pressure (enthalpy held constant). 

Parametric studies demonstrate that the pump operating characteristics can signi- 
ficantly affect the system damping. And further, it is shown that operation on positive 
pump-map slopes reduces the overall system damping. 



INTRODUCTION 


A nuclear -rocket engine in space probably will use no special startup equipment. 
Instead, it will rely on a bootstrap starting procedure. The bootstrap technique has been 
used successfully to start chemical-rocket engines, but there are significant differences 
in the nuclear -engine bootstrap. 

One difference is the much longer time for the nuclear rocket to be brought to power. 
Also, only the hydrogen propellant is used, and it follows a more complex flow path in 
the nuclear rocket. In addition there is a lack of precise knowledge of two -phase hydro- 
gen flow and boiling phenomena. All these factors caused concern that intolerable flow 
oscillations might occur during the nuclear -engine startup. Therefore, full-scale exper- 
imental tests of the bootstrap of a nuclear -rocket engine were undertaken. 

In this report, an analog computer model is used to study the mechanisms which gov- 
ern the low-frequency flow dynamics observed with liquid hydrogen flowing in two different 
nuclear -rocket cold-flow engine systems (the nuclear-rocket engine with an axial pump 
and with a radial pump) . A previous analytical simulation of the entire engine system 
with an axial pump resulted in a good overall representation (ref. 1). The model, how- 
ever, was unable to simulate the observed low-frequency flow dynamics. It was concluded 
in reference 1 that this was due mainly to the lack of an adequate two -phase equation of 
state. The purposes of this study, therefore, are (1) to generate an adequate two-phase 
state equation of the form P = P (p, h), and (2) to formulate a modeling technique which 
predicts the low-frequency flow dynamics observed in the experimental tests. 

The goal of previous analytical treatments has been to develop a mathematical model 
to simulate the nuclear -rocket engine startup. The system was studied by dividing the 
model into fixed spatial lumps, each representing a specific hardware component of the 
research apparatus shown in figure 1. These lumps were then grouped (as in fig. 2) into 
a system block diagram composed of a pump, feedline, nozzle-coolant tubes, reflector, 
core, thrust nozzle, turbine -bleed line, turbine -power-control valve, and turbine. The 
lumping was done in this manner to take advantage of system geometry and data sensor 
locations and also to limit the amount of necessary analog computing equipment. 

The differential equations describing each lump were written and then programmed 
on the analog computer. In general, the experimental data were matched to within 10 to 
25 percent except during the periods of large flow transients (ref. 1). This quasi- 
steady-state model was unable to simulate the observed low-frequency flow transients. 
There was a lack of a convenient expression relating pressure to a function of density 
and enthalpy for the two-phase fluid. This forced the dropping of certain time deriva- 
tives in the two-phase lump equations. The low-frequency dynamics that are character- 
istic of the two -phase lump were thus removed. 
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Figure L - Axial-pump nuclear -rocket cold-flow test facility. 



Figure 2. - System block diagram. 


In this study of the low-frequency flow dynamics the analog computer model used 
is a special simplication of the model used in reference 1. The present model is charac- 
terized by three fluid-lumps: a liquid lump, a two -phase lump, and a gas lump. The 
three fluid-lumps represent the system from the pump discharge to the core exit. 

The concept of the model is based on the observed state of the fluid throughout the 
system during the course of an experimental run. With liquid hydrogen flowing into the 
system, the material of the components begin to cool by giving up latent heat to the fluid. 
Throughout every experimental test, a definite, but changing, portion of the system con- 
tained two -phase fluid. The preceding portion was all liquid, while the remaining portion 
downstream was all gas (see fig. 3) . As the run proceeded and the system cooled, the 
two-phase section moved downstream. 

Differential equations are written for a liquid fluid-lump, a two -phase fluid-lump 
and a gas fluid-lump. In contrast to a conventional lump fixed in space, each fluid-lump 
is characterized by one or two variable boundaries which may shift through space (the 
system) . The positions of these moving boundaries are calculated in high-speed, repet - 



Figure 3. - Simplified model showing variation of fluid state throughout system. 


4 

















itive operation on the analog computer. 

The three -fluid-lump system is bounded on the upstream side by a run tank, inlet 
feedline, and pump. The speed of the pump is a programmed time function based on 
experimental data. The downstream side is bounded by a choked thrust nozzle. The flow 
through the turbine-bleed line is assumed negligible compared with the choked-thrust- 
nozzle flow. A block diagram of the simplified system is shown in figure 4. 

Two engine systems were studied, one with an axial pump, and the other with a radial 
pump. For each system a controlled experimental test was conducted in which step-like 
disturbances were programmed in pump speed. The resulting response in pump-inlet 



Figure 4. - Simplified system. 


flow is shown for the axial-pump system in figure 5(a) and for the radial-pump system 
in figure 5(b). For the axial-pump system the speed steps occurred every 2.0 seconds, 
with the initial step at about t = 16.0 seconds. The initial speed step for the radial- 
pump system occurred at about t = 10.5 seconds, with the subsequent steps spaced every 
3. 0 seconds. 

In general, the initial flow responses are the most oscillatory. And as time in- 
creases, the flow responses become more damped. It will be shown that this transition to 
a more damped response is due to three major factors: (1) the pump operating point, (2) 
the amount of two -phase fluid in the system, and (3) the general cooldown of the system 
material, with the subsequent lower temperature gas powering the turbine. 

To study this phenomenon, three specific flow responses are simulated. Two are for 
the axial -pump system and one is for the radial -pump system. The axial -pump system 
simulations, are of (1) the most oscillatory response, from t = 15.9 to 18.2 seconds, and 
(2) the more damped response from t = 22.2 to 24.2 seconds. Thus, the transition to a 
more damped response can be studied. Then by simulating the radial-pump system re- 
sponse (fig. 5(b)) from t = 13. 5 to 15.8 seconds, the generality of this phenomenon be- 
tween engine systems is established. 

There are two flow responses not simulated for the axial-pump system between 
t = 18.2 and t = 22.2 seconds, as illustrated in figure 5(a). During these responses the 
pump operated both in-stall and out-of-stall. There is some uncertainty about the mech- 
anisms involved during the transition into and out of stall. It was decided not to introduce 
this uncertainty into the simulation. Hence, the first time period during which the pump 
operated completely out-of-stall (t = 22.2 to 24.2 sec) was selected. 
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Figure 5. - Pump-inlet flow transients resulting from controlled, step-like, pump speed disturbance 
in experimental test. 


6 



During the three simulated time periods, the two -phase fluid was located in approxi- 
mately the same section of the system. The saturated-liquid fluid boundary was inside the 
nozzle -coolant tubes and the saturated-gas fluid boundary was just inside the reflector. 

This choice of time periods, therefore, also simplified the computer model. The high- 
speed repetitive-operation portion of the analog, which locates the moving fluid boundaries, 
need only consider the system from the nozzle -inlet manifold to the reflector exit. The 
result is a large reduction in necessary analog computing equipment. 

In the remainder of the report the two nuclear -rocket cold-flow engine systems will 
be briefly described. The analog computer model will be presented, with a discussion on 
mass-flow storage and with most equation derivations referred to appendixes A to E. The 
analytical data from the system response simulations of the three time periods of interest 
will then be discussed and compared with experimental values. 


DESCRIPTION OF NUCLEAR-ROCKET ENGINE SYSTEMS 
The Axial-Pump System 

The nuclear -rocket cold-flow test facility was located at the Plum Brook Station of 
the Lewis Research Center. The altitude exhaust system (fig. 1) maintained a rocket- 
nozzle -exhaust pressure of approximately 0. 5 psia (3450 N/m ). The rocket engine was 
mounted in the structure at the top right in the figure in a down-firing position. 

A schematic diagram of the rocket-engine test package using the axial pump is also 
shown in figure 1. The test hardware consisted of a turbopump assembly, unfueled 
reactor, and a supersonic exhaust nozzle. The liquid-hydrogen run tank which fed the 
pump, had a servocontrolled pressurizing system. An 8-inch (0. 2-m) diameter line, 
which was equipped with a turbine flowmeter, connected the pump to the run tank. 

The turbopump assembly consisted of a Rocketdyne MARK -IX turbopump, turbine- 
power-control valve, and turbine-bleed line The liquid-hydrogen pump was composed 
of an axial-entrance mixed-flow axial -discharge inducer stage, six identical high- 
pressure axial -flow stages, and a single-outlet collecting volute. The turbine was a six- 
stage pressure -compounded axial-flow unit designed for operation with hot or ambient 
temperature gaseous hydrogen. 

The pump discharged into a 4-inch (0. 1-m) diameter stainless steel flight -weight 
feedline equipped with a servocontrolled butterfly valve used for flow control. Approx- 
imately 18 feet (5.5 m) downstream from the main flow-control valve, the propellant 
line divided equally into three smaller ducts that fed the nozzle -inlet manifold. 

The regeneratively cooled (with liquid hydrogen) tubular -wall nozzle utilized single - 
pass cooling. The liquid hydrogen entered the tubes through three equally spaced inlet 
connections on the manifold and flowed toward the reactor end. 
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On discharging from the nozzle -coolant tubes, the fluid entered the reflector-inlet- 
manifold section of the reactor assembly. From there it passed through the reflector 
component that was composed of an inner graphite reflector, an outer aluminum reflec- 
tor, simulated aluminum control rods, and the aluminum pressure vessel. The fluid 
then entered the reactor core, flowing through the many parallel coolant passages in the 
unloaded graphite reactor core. 

The warm hydrogen emerged from the core into the nozzle thrust chamber where 
most of it was exhausted into the facility vacuum exhaust. A small portion of the fluid, 
however, was bled to the turbine through a 3-inch (0.076-m) diameter bleed port, the 
turbine-bleed line, and the servocontrolled turbine-power-control valve. Altogether, 
the rocket-engine system tested was a full-scale model of what might have been a flight 
version of a system using the Kiwi reactor. 


The Radial-Pump System 

This nuclear -rocket cold-flow test facility was also located at the Plum Brook 
Station. Except for a few specific changes, this test facility (ref. 2) was assembled 
identically as the test facility for the axial -pump system (fig. 1) . 

The turbopump was changed to an Aerojet Mark m, Mod 4 turbopump. The liquid- 
hydrogen pump was a radial flow, centrifugal type with an integral three-blade helical 
inducer. The turbine was a two-stage pres sure -compounded impulse type. 

The only other change was in the propellant ducting. The inlet piping that connected 
the pump to the run tank was lengthened by about 6 feet (1.8 m) . 


AN ANALOG COMPUTER PROGRAM FOR SIMPLIFIED ENGINE SYSTEMS 

An analog computer program is developed for the nuclear-rocket cold-flow engine 
systems described in the previous section. The purpose of this program is to study 
system transient response to step-like, pump-speed disturbances. Experimental data 
are utilized in several areas of the analog model to simplify the program and help to 
establish the operating levels of pressure and temperature observed in the experimental 
tests. Emphasis is placed on analyzing the frequency and damping content of the computer 
results, rather than on how well the experimental data are reproduced. 
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Principles of Operation 


The simplified representation for the systems simulation is shown in figure 6. 
Experimental data are utilized for the pump speed and for the nozzle -inlet-manifold fluid 
enthalpy. The pump-speed forcing function, combined with neglecting turbine-bleed 
line flow, eliminates the need to simulate the turbine -bleed line, servocontrolled 
turbine -power -control valve, and turbine. The enthalpy forcing function defines a ref- 
erence point in a section of the analog program called the spatial-distribution model. 
This, as demonstrated in appendix E, further simplifies the computer representation. 

The analog computer program consists of two parts: (1) a slow-speed analog, and 
(2) a high-speed analog. With the exception of the spatial -distribution model, the entire 
system in figure 6 is programmed on the slow-speed analog. The spatial-distribution 
model is solved on the high-speed analog. 



Figure 6. - Simplified system for analog computer program. 


The slow-speed analog is operated at 1/100 real time; that is, 100 seconds of com- 
puter operation are required to generate one second of solution. The high-speed analog 
is operated repetitively at 5 hertz relative to slow-speed analog operation. It follows 
from these two rates that the spatial-distribution model updates its outputs 500 times per 
second of solution. 

The high-speed analog solution is assumed to take place fast enough such that the 
input variables from the slow-speed analog appear as constants. Spatial distributions 
are computed with the aid of high-speed integration with distance x as the independent 
variable. Certain variables (see appendix D) are then sampled, stored, and then sent 
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as updated inputs to the slow-speed analog. The high-speed analog is then reset and the 
computation restarted with new inputs from the slow-speed analog. 


Model Lumping 

A diagram of the model lumping was presented in figure 6. The first lump to be 
considered is the run tank. For each engine system, the run tank is treated as a 
constant -pressure liquid-hydrogen supply. 

The all-liquid pump-inlet line is a stationary, constant-geometry lump. It is as- 
sumed that the walls of the inlet line have lost all of their latent heat content during a 
pretest cooldown. Hence, the operation of this lump is based on the conservation of fluid 
mass and momentum principles. Fluid energy considerations are neglected. 

A normalized, steady-state pump characteristic is used to represent the pump. A 
comparison of the operating characteristics of the two types of pumps is shown in figure 7. 



Figure 7. - Operating characteristics of pumps used in axial- 
and radial-pump systems. 
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The data for the radial-pump operating characteristic was obtained from a series of 
special pump tests. However, the axial-pump operating characteristic was derived with 
quasi-steady-state data points from the series of tests conducted. 

Consider the portion of the system in figure 6 between the pump and the thrust 
nozzle. The mathematical model for this section consists of three parts: (1) the three- 
fluid-lump model, (2) the spatial-distribution model, and (3) the system-heat-transfer 
model. Arrows are drawn in the figure to illustrate the general flow of information 
between the three parts. The fundamental principle in the operation of this model con- 
cerns the storage of mass flow in the system. It is shown in appendix A that the amount 
of storage is proportional to the thermodynamic fluid property [3p/9P] ^ and, further, 
that most of the storage takes place in the two-phase fluid near the saturated-liquid 
boundary. 

The three-fluid-lump model consists of three variable, special lumps called fluid- 
lumps: a liquid lump, a two-phase lump, and a gas lump. The two-phase fluid is treated 
as a homogeneous mixture; the state equation used is derived in appendix B. The opera- 
tion of this three-fluid-lump model is based on the conservation of fluid mass and momen- 
tum principles. 

The primary function of the spatial-distribution model is to maintain the conservation 
of fluid energy . The operation of this lump is characterized by high-speed, repetitive 
generation of the x-distributions of certain variables. The x-distribution of the fluid 
property [3p/9 P]^, for instance, is calculated and averaged every cycle. This, in turn, 
provides for approximating the amount of mass-flow storage in a fluid-lump (appendix A) . 

The system-heat-transfer model consists of several stationary lumps, each repre- 
senting a fixed amount of system mass. The operation of this part is based on the con- 
servation of energy as applied to the system material. Asa result, approximate cool- 
down rates are determined for different sections of the system. 

The combined operation of the three-fluid-lump, the spatial-distribution, and the 
system -heat-transfer models fully represents the section of the engine systems between 
the pump and the thrust nozzle. In summary, the main fluid-dynamic calculations are 
performed by the three -fluid-lump model, with the system-heat -transfer model supplying 
the more slowly varying heat-transfer dynamics. The spatial-distribution model 
performs high-speed averaging and x-distribution calculations. 

The thrust nozzle represents the final lump (fig. 6) to be considered. The nozzle is 
assumed choked at all times. The thrust-nozzle -chamber pressure calculation is based 
on the conservation of fluid mass between the core exit and the nozzle exhaust. The 
turbine-bleed-line flow is neglected. Neglecting energy considerations, the core-exit 
fluid temperature is used in the thrust-nozzle chamber. 

The basic mathematical equations programmed on the analog computer are derived 
in appendixes A, B, and C. Specific equations are developed for each lump in appen- 
dixes D and E . 
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DISCUSSION OF RESULTS 


The analytical data presented in this section are the results of dynamics studies using 
a simplified model for the engine systems. Experimental data were programmed into 
those areas which would add unnecessary complexity to the model and have a negligible 
effect on the low-frequency dynamics under study. The emphasis of the following dis- 
cussion, therefore, shall be mainly on the frequency and damping content of the computer 
data, rather than on how well the experimental data are reproduced. 


Input Forcing Functions 

Pump-speed forcing functions. - The pump-speed forcing functions used are illus- 
trated in figure 8. The curves show the response of the pump to step-function speed 
demands initiated in the turbine -power -control valve control system. The control-system 
function was to rapidly accelerate or decelerate the pump speed between 4500 and 5500 rpm 
to the demanded speed. 

The speed curve in figure 8(a) shows the response in the axial-pump system to a 
down-speed step demand. The demand, initiated at t = 16. 3 seconds, called for a decel- 
eration from 5500 to 4500 rpm. The preceding 0.4 second is included in the forcing 
function, however, to establish a low flow acceleration starting point. To minimize 
computation starting transients, it is desirable to begin a computer solution at a steady-, 
or at least a qua si -steady-, state point. During the 0.4 second from t = 15.9 to 16.3 
seconds the turbine-power-control valve control system was decelerating the pump speed 
toward 5500 rpm to compensate for the approximately 350-rpm overshoot existing at 
t = 15.9 seconds. 

The speed curve for the second time period considered in the axial-pump system is 
presented in figure 8(b). The speed demand initiated at t = 22.2 seconds called for an 
acceleration from 4500 to 5500 rpm. At this time in the test, however, the lower-energy 
turbine-inlet gas was unable to accelerate the turbine sufficiently to satisfy the demand. 
Consequently, as shown in figure 8(b), the speed response fell short to the demanded 
5500 rpm and, further, the turbine -power -control system was unable to maintain a steady 
final value of pump speed. This more sluggish system response has an appreciable effect 
on the subsequent flow dynamics, as will become evident later. 

The speed curve in figure 8(c) shows the response in the radial-pump system to a 
step demand. The demand, calling for a deceleration from about 5500 to 4500 rpm, was 
initiated at t = 13.5 seconds. The resulting speed response was similar to the axial- 
pump system response shown in figure 8(a). The radial -pump system response, however, 
has a small oscillation superimposed on the general deceleration waveform. This more 
oscillatory system response is also apparent in the resulting low-frequency flow dynamics, 
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as will be demonstrated later. 

Nozzle-inlet-manifold fluid - enthalpy forcing fu nctions. - The curves shown in fig- 
ure 9 are the forcing functions used for the nozzle-inlet-manifold fluid enthalpy h n - . The 
experimental data points were calculated from measured values of fluid temperature and 
pressure. The liquid-hydrogen -properties digital subroutine (presented in ref. 3) was 
used. 

The solid line in each figure represents a linear least-squares curve fit of the exper- 
imental data points. The dotted line is an approximate saturated-liquid-enthalpy refer- 
ence line and is determined by the nozzle-inlet-manifold pressure. The vertical distance 
that the solid line is below the dotted curve indicates the degree of fluid subcooling. Two- 
phase fluid exists whenever the solid line is above the dashed reference curve. 

The use of a forcing function establishes a reference point in the system near 
the saturated-liquid fluid boundary. This eliminates the need for fluid and material 
energy equations to predict the fluid enthalpy spatial distribution for the upstream, low- 
heat-content portion of the systems. The result is a significant reduction in the amount 
of necessary analog computing equipment. 

The enthalpy curves in figure 9(a) are for the axial -pump time period beginning at 
t = 15.9 seconds. Two-phase fluid exists in the nozzle-inlet manifold from t = 16.8 to 
17.7 seconds. No experimental data points are presented for these times since fluid 
temperature and pressure do not define a unique value of two-phase fluid enthalpy. Rela- 
tive to the curves in figures 9(b) and (c), the fluid enthalpy in this figure shows the least 
amount of subcooling. As will be shown later, the amount of fluid subcooling in the 
nozzle-inlet manifold significantly affects the low-frequency flow dynamics, since it is a 
major factor in determining the amount of two -phase fluid in the system. 

The largest amount of fluid subcooling is indicated by the enthalpy curves in fig- 
ure 9(b). These curves are also for the axial-pump system, but for the later time period 
beginning at t = 22.2 seconds. At this point in the test, the rate of heat addition to the 
hydrogen from the cooled, upstream system material was low. This is illustrated by the 
small time rate of change of the h ni curve. 

The enthalpy curves in figure 9(c) are for the radial-pump time period beginning at 
t = 13. 5 seconds. The significant amount of fluid subcooling is similar to that shown in 
figure 9(b). As indicated by the relatively large time rate of change of the h ni curve, 
the upstream system material was rapidly cooling off during this period of the test. 
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(c) For radial -pump system time period. 

Figure 9. - Programmed nozzle-inlet-manifold fluid enthalpy forcing function. 
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Comparison of Predicted and Measured Data 


In general, the predicted fluid flow, pressure, and temperature data are within 
±13 percent of the measured values, as illustrated by figures 10 to 14. The predicted 
values are identified in the figures with solid lines, and the experimental measurements 
are identified with symbols. 

The predicted and measured values for pump-inlet flow are compared in figure 10. 

In figure 11, the predicted and calculated data for the thrust-nozzle-chamber flow are 
compared. The experimental thrust-nozzle -chamber flow is calculated with a choked- 
nozzle equation (appendix E) from experimental measurements of thrust -nozzle -chamber 
fluid temperature and pressure. The pump-inlet and thrust -nozzle -chamber flows togeth- 
er show that the proper mass-flow storage in the system has been approximated. 

The frequency and damping content of the predicted pump -inlet flow agree reasonably 
well with observation. The predicted period of oscillation for the solution beginning at 
t = 15.9 seconds in figure 10(a) shows good agreement, except for a slight time shift of 
approximately 0. 1 second. The largest error in the flow period is shown in figure 10(b) 
during the initial portion of the transient. Discounting the approximate 0. 1-second time 
shift, the period is off by about 13 percent at t = 14. 3 seconds. 

Figure 12 shows the pump-inlet and pump-outlet pressures. The predicted values of 
nozzle -inlet -manifold pressure and the thrust-nozzle -chamber pressure are compared 
with the experimental values in figure 13. In general, the calculated pressure levels 
compare well with observation. The pressure rise between the pump outlet and inlet 
(fig. 12) is governed by the pump operating characteristics. The nozzle-inlet-manifold 
pressure (fig. 13) is an initial condition for the spatial -distribution model solutions and 
is approximated with a simplified equation (appendix D) . 

The predicted thrust -nozzle -chamber fluid temperature is compared with an experi- 
mental average value in figure 14. For the axial-pump system the experimental temper- 
ature represents an arithmetic average of the five available measured values. The ex- 
perimental temperature for the radial-pump system, however, involves additional meas- 
ured values, plus a weighting technique to account for the cross-sectional distribution of 
temperature. In general, the results are close to observation, considering the single- 
lump treatment of the entire gas phase . 

The two metal temperatures generated in the computer model are shown in figures 15 
and 16. The temperature T w ng is an average value for the nozzle -coolant-tubes shell 
material and is presented in figure 15. The indicated experimental values are measured 
about midway along the shell length. As shown, the predicted cooldown of this material 
agrees very well with experimental data. 

The metal temperature T w ri _ ce (fig- 16) is an average value for the entire gas 
fluid-lump. The experimental values shown were measured at six designated points be- 
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(b) For radial -pump system time period. 
Figure 10. - Pump-inlet flow. 
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(b) For radial-pump system time period. 
Figure 11. - Thrust-nozzle-chamber flow. 
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Figure 12. - Pump-inlet and pump-outlet pressures. 
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(a) For axial-pump system time period beginning at t* 15.9 seconds. 
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Figure 13. - Nozzle-inlet-manifold and thrust-nozzle-chamber pressures. 
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Figure 14. - Thrust-nozzle-chamber fluid temperature. 
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Figure 15. - Nozzle-coolant -tubes shell average metal temperature. 
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Fioure 16. - Reflector-inlet-to-core-exit average metal temperature. 
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tween the reflector inlet and the core exit. In general, the predicted values are low. 
However, the slopes of the curves indicate that the observed cooldown rates have been 
approximately duplicated. Hence the total input heat flux to the gas fluid-lump has 
been approximated. The value of Q g , in turn, is used to calculate the core-exit fluid 
temperature (fig. 14). The equations (pp. 62 and 73, appendix D), which demonstrate this 
effect are as follows : 


3T. 


w, ri-ce _ 

at 


Q 


g 


M ri-ce c p, sg-ce 


Q s 

T, = T , + s 

b, ce sat ^ c 


g P,g 


(Symbols are defined in appendix F . ) 

Predicted Movements of Fluid Boundaries 

The curves of figure 17 represent the calculated movements of the saturated-liquid 
and the saturated-gas fluid boundaries. The values plotted correspond to the respective 
distances of the fluid boundaries measured downstream from the nozzle-inlet manifold. 
No experimental measurements, unfortunately, are available for comparison. However, 
a few qualitative comments may be made. 

In general, the saturated-gas fluid boundary moves less than the saturated-liquid 
fluid boundary. This downstream fluid boundary is situated inside the high-input - heat- 
flux reflector throughout all of the computer solutions. Thus, only a small spatial shift 
of the boundary is required to compensate for an impinging flow disturbance. The effect 
is exemplified in the axial -pump system simulation by comparing the boundary move- 
ments computed during the two speed disturbances (figs. 17(a) and (b)). The largest 
movement is observed during the second disturbance, after the reflector material has 
cooled to a lower temperature. 

The nozzle-coolant-tubes region upstream of the throat is a low-input - heat-flux 
region, and the remaining portion of tubes downstream is a relatively high-input - heat- 
flux region, as illustrated in figure 18. Accordingly, the most significant spatial move- 
ments of the saturated-liquid fluid boundary takes place upstream of the throat, as 
shown in figure 17. The time rate of change of the boundary position clearly decreases 
(fig. 17(b)) as the boundary shifts downstream of the throat. 


24 



Distance, 


2. i i — 7 r 


Saturated gas 


Reflector inlet 


Saturated gas 


1.5k ° 5 


Reflector inlet 


Saturated liquid 


.Saturated liquid 


ol I N. 1 I I ^^1 , 6 L zl I I I I I I 

15.8 16.2 16.6 17.0 17.4 17.8 18.2 22.0 22.4 22.8 23.2 23.6 24.0 24.4 

Time, t, sec Time, t, sec Time, t , sec 

(a) For axial-pump system time period beginning at t c 15. 9 seconds. (b) For axial-pump system period beginning at t = 22. 2 seconds. 

2.1 1 — 7 r~ Saturated gas 


Reflector inlet 


•6 2h 



.31 H I I I 1 _J | 

13.4 13.8 14.2 14.6 15.0 15.4 15.8 

Time, t, sec 

tc) For radial-pump system time period. 

Figure 17. - Predicted distance downstream of nozzle-inlet-manifold of two-phase saturated fluid boundaries. 
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Figure 18. - Simplified diagram of nozzle -in let -manifold -to -ref lector -exit region 
showing areas of high and low wall temperatures. 

Effect of the Pump Characteristics 

Operation on positively sloped segments of the normalized pump characteristics 
(see fig. 7) will reduce the damping of flow transients. In fact, it is concluded in ref- 
erence 4 that overall system instability will occur if the pump -map slope exceeds a 
certain critical, positive value. To demonstrate the sensitivity of the pump -map slope 
and the tendancy of positive slopes to counteract friction losses, parametric studies for 
the two engine systems will be presented. The resulting curves, therefore, are para- 
metric solutions and not actual system transients. 

Axial-pump sys tem. - The pump operates in stall (positive slopes) during the entire 
solution beginning at t = 15.9 seconds and operates out-of-stall (negative slopes) during 
the solution beginning at t = 22.2 seconds. An approximate steady-state pump map 
was derived for the out-of-stall region using data from 10 experimental tests. No steady- 
state data points were found for the stall region, however. Consequently, normal opera- 
tion in-stall is assumed to take place on a constant, positively sloped characteristic. A 
value for the stall slope was selected which yielded solutions in good agreement with 
experimental data. 

Six constant -slope variations of the pump characteristic are assumed. The three 
in-stall variations and the three out-of-stall variations are shown in figure 19. A portion 
of the negatively sloped, out-of-stall pump map is included also as a dotted line. For 
each set of three variations considered, the solution is begun at the same initial conditions, 
and the subsequent pump operation is defined by the constant -slope variations. 

The curves in figure 20 represent the solutions for pump-inlet flow, which result 
from the parametric variation of the pump -map slope. The flows resulting from the in- 
stall variations are shown in figure 20(a), and the out-of-stall curves are shown in fig- 
ure 20(b). Comparing the in-stall and out-of-stall solutions thus establishes the high 
sensitivity of the system response to a positive pump-map slope. And further, the de- 
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Figure 19. - In-stall and out-of-stall pump-map 
slope variations considered for axial-pump 
system time periods. 
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crease in solution damping supports the theory that overall system instability can result 
if the slope of inlet -pumping characteristics becomes too positive (ref. 4). 

Radia l -pump system . - The pump operates out-of-stall during the speed disturbance 
considered. The steady-state pump map was derived from a series of special pump tests, 
and a portion of the out-of-stall region is compared with the axial pump map in figure 21. 
As illustrated in the figure, a particular characteristic of the radial pump is the presence 
of positive slopes in the out-of-stall region. The solution (beginning at t = 13. 5 sec) 
operates in the positively sloped region of the map. 

The sensitivity of the system to the slope of the pump map is similar to that observed 
for the axial -pump system. To demonstrate this, a pump -inlet-flow solution using a zero 
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Figure 21. - Positive-slope pump-map variations for axial- 
and radial-pump systems. 
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Time, t, sec 


Figure 22 . - Radial-pump system pump-inlet flow solutions using 
zero and regular pump-map slopes. 


pump-map slope is compared in figure 22 with a solution using the regular pump map. 

The regular -pump -map solution is more oscillatory due to the positive slopes of the pump 
map. 

The difference between the regular -pump -map solution and the zero-pump-map so- 
lution, however, is much smaller for this engine system than it was for the in-stall 
axial -pump system (fig. 20(a)). The difference is mainly attributed to the more oscilla- 
tory pump -speed forcing function which is used with the radial pump. As presented 
earlier while discussing figure 8, the radial-pump system response to the speed demand 
was less damped. The resulting pump -speed waveform, shown in figure 8(c), has a 
small oscillation superimposed on the general smooth trend. 

Another factor which contributes to the more oscillatory radial-pump system re- 
sponse is the longer pump-inlet feedline. The effect of the longer line is to increase 
the amount of fluid (liquid) inertia. However, this reduced-damping effect is small 
because of the large diameter of the pump-inlet feedline and the relatively large amount 
of liquid inertia existing downstream of the pump. 
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Importance of Correct Averaging 


When the modeling technique for a system uses gross spatial lumps as approxi- 
mations to the spatial dependence of the describing partial differential equations, the 
type of lump averaging becomes important. In particular, when treating the mass-flow 
storage in the nuclear -rocket cold-flow engine systems described, a special type of 
averaging technique proved necessary. 

It was found in the analysis in appendix A, that to approximate the mass-flow stor- 
age in a gross fluid-lump it is necessary to calculate the average value of the thermo- 
dynamic fluid property [0p/0P]^. Simulation of the observed low-frequency flow dy- 
namics, in fact, hinges also on the continued updating of this average. Hence, the 
average values for pressure, density, and enthalpy in the two-phase fluid-lump are 
related as follows : 


and 


P 2p P ^2p,av’ h 2p, av^ 


h 2p, av 




(Throughout the equations in this report, parentheses are used solely for functional 
notation, and brackets and braces are used for all other purposes.) A set of curves 
for the second expression is derivable from the hydrogen-properties tables in refer- 
ence 3. 

To demonstrate the importance of calculating the correct average for h 2 p ay , two 
basic averaging schemes are compared: (1) relating h 2 p av to the arithmetic spatial 
average of [3p/9P] and (2) calculating the arithmetic spatial average of h 2 p ay . For 
each scheme, the effect of holding the average at its initial value during the solution is 
also considered. The solutions shall be referred to as those resulting from the following 
enthalpy -averaging techniques: 

(1) The average two-phase fluid enthalpy h 2p ay predicted by the previous equa- 
tions, with the average thermodynamic fluid property updated continually in high-speed, 
repetitive operation by the equation 
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calculated once and then 


(2) The hgp ay predicted as in (1) but with [9p/3P] h ay 
held constant throughtout the computer solution. 

(3) The hgp ay predicted by the equation 


/ 

•/v 


Sg 


2p,av 


_ X sZ 


h 2 p(x)dx 


sg si 


continually in high-speed, repetitive operation. 

(4) The hgp ay predicted as in (3) but with the average held constant at its initial 
value throughout the computer solution. 

The results for both engine systems are consistent. The solution for pump-inlet 
flow which is closest to experimental data uses enthalpy-averaging technique (1), that is, 
the technique which helps establish the correct mass -flow storage in the system. The 
curves in figure 23 show the results of using the four averaging techniques for the down- 
speed-step radial-pump system simulation. The results for the axial-pump system are 
similar. 

Comparing figure 23(c) with 23(a) demonstrates the effect of using an arithmetic 
spatial averaging technique for h 2p,av‘ The curve in figure 23(c) is characteristic 
of a system with less capacitance, or storing capability. The period and the amplitude 
of oscillation are both reduced. 

The importance of continually updating hgp ay is evident by comparing figure 23(a) 
with 23(b) and figure 23(c) with 23(d). The solution in figure 23(b) is continuous; the 
plotting was merely discontinued and then resumed at the x-axis. The hg p ay from 
technique (1) never changed more than 16 percent from its initial value during the total 
solution. 

The predicted mass-flow storage in the system is obtained by taking the difference 
between the predicted values of pump-inlet flow and thrust-nozzle -chamber flow. The 
curves in figure 24 show the predicted mass -flow storage which results from using each 
of the four enthalpy-averaging techniques. The calculated values shown are determined 
from experimental values of pump-inlet flow and calculated values of thrust-nozzle- 
chamber flow. These results are presented, for the positive-speed-step axial-pump 
system simulation, since they most clearly illustrate the effect of the enthalpy -averaging 
technique on the system mass-flow storage. The best solution (shown in fig. 24(a)) uses 
enthalpy-averaging technique (1) . Whereas an obvious lack of mass-flow storage is 
illustrated in figure 24(c) . 

By comparing figure 24(a) with 24(b) and figure 24(c) with 24(d) the significance of 
continually updating the average two-phase fluid enthalpy becomes obvious. During the 
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solution, the hgp av from technique (1) never changed more than 3 percent from its 
initial value. 


Additional Effects 

Nozzle -inlet-manifold fluid enthalpy . - An important forcing function is the nozzle- 
inlet-manifold fluid enthalpy h n ^. The magnitude of h R . has a large effect on the com- 
puter solution since its value defines a reference point for the fluid-enthalpy x- 

distribution calculations . Consider the errors AT, „ and AP„ in the measured values 

b, e e 

of nozzle -inlet-manifold fluid temperature and pressure, respectively. The resulting 
error in the fluid enthalpy would be 


Ah e = 


_3h 

0T 


AT 


b, e + 


_0h 

[_0Pjrp 


AP, 


For liquid-hydrogen fluid, the more dominant error term is that due to the fluid tem- 
perature error, such that 


Ah e 



b, e 


The curves in figure 25 show the appreciable effect on the pump -inlet-flow solution 
of 0. 5 and -0. 5 Rankine degrees (±0.278 K) errors in the measured value of the nozzle- 
inlet-manifold fluid temperature. The higher temperature fluid results with a value of 
h n i closer to the value of saturated-liquid enthalpy. Hence, a higher value of T^ 
necessarily corresponds to a larger two -phase volume because of the relatively low 
input heat flux of the nozzle-coolant tubes. 

It is shown in reference 4 that the presence of a boiling, two-phase compressible 
volume can reduce the system stability and make it more oscillatory. The curves in 
figure 25 support this, since the most oscillatory solution is obtained with the warmest 
fluid. The relatively large effect of the temperature change, however, was not expected. 

Use of steady -state pump map . - The pump operating point is assumed defined at all 
times on a steady-state pump map. In general, this may not be the case, since the solu- 
tions contain large transients. Significant deviations from the steady-state characteris- 
tics can take place under certain transient conditions. Based on the response sensitivity 
observed in the pump -map-slope parametric variations, the introduction of pump dynam- 
ics into the model could have a meaningful effect of the solutions. 
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SUMMARY OF RESULTS 


Simplified models of two nuclear-rocket cold-flow engine systems were programmed 
on an analog computer to study the mechanisms governing the observed low-frequency 
flow dynamics. A summary of the results of this study are as follows: 

1. An adequate two-phase state equation for the case when thermodynamic equilib- 
rium may be assumed between the gas and liquid phases is 

P 2p * p 2p[“< P 2p> + « W 


where P 2p is the two -phase fluid pressure, p 2p is the two -phase fluid density, h 2p is 
the two -phase fluid enthalpy, and a and /3 are coefficients which are functions of 
pressure. 

2. The key to obtaining a low-frequency response is to obtain the proper mass-flow 
storage. This is achieved by relating the average fluid enthalpy of a fluid-lump to the 
average value of the thermodynamic property [dp/d Pj^; that is. 


av 




P av P ^av’ h av^ 


where [ 3p/3P] ^ is the partial derivative of density with respect to pressure (enthalpy 
held constant), P is the fluid pressure, p is the fluid density, h is the fluid enthalpy, 
and the subscript av indicates an average quantity. 

3. Analysis shows that most of the mass-flow storage occurs in the region of low 
mass quality, that is, in the two -phase fluid in the neighborhood of the saturated-liquid 
boundary. 

4. The pump operating-characteristics can significantly affect the system damping. 
In fact, operation on positive pump-map slopes reduces the overall system damping. 

The following remarks are made concerning the analytical model used: The mag- 
nitude of nozzle-inlet-manifold fluid enthalpy is important since it defines a reference 
point for the fluid-enthalpy x -distribution calculations. Introducing the possible exper- 
imental error in determining this forcing function results in significant changes in the 
solution for pump -inlet flow. 

A small amount of computational instability occurred with the combined operation of 
the slow-speed analog and the high-speed analog. The output from the high-speed track- 
and-store circuits was noisy and had to be filtered. It was found that the inherent errors 
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associated with the track-and-store units (±0. 5 percent) would cause a sampled variable 
to become noisy. Whenever such an unfiltered variable was used in a relatively high- 
gain (X100) circuit, the entire simulation was unstable. The size of the required noise - 
filter capacitors, of course, depends on the repetitive operation rate. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 11, 1969, 

126-15-02-36-22. 
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APPENDIX A 


DERIVATION OF VARIABLE-VOLUME EQUATIONS 

Consider a gross, one -dimensional system lump, bounded by inlet and exit surfaces 
(which may be variable) and composed of l equal sublumps at any instant of time (see 
fig. 26). It is desired to approximate the proper mass-flow storage, pressure drop. 


12 3 4 


' 

. . . 

l-l 

l 


Figure 26. - Sublumping of gross, one-dimensional system lump. 


and heat transfer for this section of the system by single equations of the respective 
forms 


aW tot=^Kv V l0t] 


AP, , = AP. + AP t + AP + AP 
tot iner fr mom gr 


Qtot - E Qi 

i=l 


Approximating Mass-Flow Storage 

Mass -flow distribution for gross lump . - Consider first the system shown in fig- 
ure 26 and assume that the bounding surfaces are fixed. The amount of mass-flow 
storage in the i n sublump is given by 


AW, = 


at 
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Hence, the total mass -flow storage would be 



3p i 


at 


The state equation of the fluid in the 


i th 


sublump may be written as 


'>i=p 1 < I V h i> 


By the chain rule. 


9p i 


at 



Rp 


LahJ 


^i 

P, i at 


The total mass -flow storage, therefore, is 



Figure 27 shows the two thermodynamic properties [dp/d P]^ and -[3p/ah]p for 100 per- 
cent parahydrogen as a function of fluid enthalpy. The curves were calculated for con- 
stant pressures of 5, 20, and 80 psia (34, 138, and 552 kN/m ) from the data on hydrogen 
properties in reference 3. The variations of the two thermodynamic properties are 
shown for the liquid, two-phase, and gaseous fluid states. The vertical dotted lines are 
located at the values of saturated-liquid and saturated-gas enthalpies. 

The pressure and enthalpy rate of change terms (the time derivatives) are assumed 
approximately uniform throughout a variable -volume lump. Therefore, we can write 
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The maximum mass-flow storage occurs at the point in the lump where [dp/d P]^ and 
[3p/3h] h are largest. Clearly, the most mass -flow accumulation will take place in the 
two-phase fluid. In addition, it is evident that most of this mass-flow storage will occur 
in regions of low mass quality, that is, in the vicinity of the saturated-liquid boundary. 

Experimental data indicate that for the system under study the enthalpy time rate of 
change term is small compared with the pressure time rate of change term. An approx- 
imate expression for the total mass -flow storage, therefore, is 


AW tot = V tot 


pp] 


3P" 

3P 

h, av 

3t 


av 




Hence, average values for pressure, density, and enthalpy will be related through the 
average value of the thermodynamic property [3p/3P]^. Mathematically, this may be 
written as 


av 



and 


^*av ^^av’ *V v ) 

Continuity equation for gross lump. - For the constant -volume approach the mass- 
flow storage may be written as 


AW, , = W. 
tot in 


- w 


out 


For the variable -volume approach, however, the inlet and exit boundaries are variable. 
If the inlet boundary moves with a velocity Ug - n and the exit boundary moves with a 
velocity Ug ou ^, the mass-flow storage equation becomes 

"tot - tw - PAU B ]. n -[W-pAU B ] out 


The continuity equation for the gross lump, therefore, may be written as 



Approximating Pressure Drop 

Consider similar sublumping as in the previous section. The pressure drop across 
the i** 1 sublump would be 


AP, . . = AP. . + AP„ . + AP . + AP 

tot, 1 iner, 1 fr, 1 mom, 1 gr, 1 


The total pressure drop across the gross lump would be 


l 

Ap tot = y> p tot.i 
i=l 

III l 

ap. .+v > ap £ . ap . +y ^ ap 

/ / iner, 1 / / fr, 1 / v mom, 1 / > gr, 1 

i=l i=l i=l i=l 

The individual terms will now be considered. Assume that each sublump has constant 
flow area (the sublumps need not be of equal length). Consider only the geometry var- 
iations between sublumps. Assume that flow is uniform. 

Total inertia pressure drop . - 

l l 

ap. = > ap. . = — ) JLJLrwAxi. 
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Total friction pressure drop. - For circular passages 
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The coefficient of friction (ref. 5) is usually defined as 


2 gr 

f = w 
p b U b 


For turbulent flow with Re < 200 000 an approximate expression is 


f _ 0. 046 
~Re 0 * 2 


= 0, 046 
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where the viscosity is to be evaluated at the film temperature and pressure. Bends and 
entrance effects may be included in an additional loss term K losg . The effective wall 
shear of the i th sublump may, therefore, be written as 
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for two -phase fluid (appendix C), where 
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“JT - + K loss, i 
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Hence, the following equation approximately represents the total friction pressure drop: 


AP fr = 
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for liquid and gas, and 
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for two-phase fluid (appendix C), where 


p av 


^ inPout 
^in + ^out 


An average value of film viscosity is assumed. 
Total memontum pressure drop. - 


AP 


mom 


AP 



[W - pAU B l 


in 


^in^in 




A more convenient and approximate form for this equation is 


AP 


mom 



It is assumed that the area change from sublump to sublump is accounted for with the loss 
term in the friction pressure-drop equation. The preceeding expression for the complete 
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momentum pressure drop, although only an approximate representation, provides a 
convenient variable coefficient for analog computer solution. 

Total gravity pressure drop. - 



n 


Momentum equation for_gross lump. - Beginning with the complete pressure drop 
equation 


AP tot AP iner + AP fr + AP mom + AP gr 


the following form of the momentum equation for the gross lump may be derived: 


aw 

at 
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tot 
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K av 
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where 


a = 


ng 
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for liquid and gas r and 



for two -phase fluid, 


c = 


' tot 


Z A i iX i 

i=l 


and 


v 


av 


1 _ ^in + ^out 

P- av ^in P out 


Approximating Heat Transfer 

Consider the same sublumping as in the previous two sections. The heat flux 
in the i** 1 sublump from the wall to the fluid is 


=*i A ht,i [T w,i - T b,i] 


The total heat flux for the gross lump is then 


«tot’ £ 


i=l 


z 

i=l 


^f.A. . . 

l ht, 1 


[ T w - T b] 
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where an average wall-to -fluid temperature difference is assumed. 

In this study heat transfer as a gross lump will only be considered in the gas fluid- 
lump. The heat-transfer coefficient is determined from the conventional correlation 
in reference 5; that is, 


Nu b = 0.023 Re b ' 8 Pr b ’ 4 

where the fluid properties are evaluated at the bulk fluid temperature. The heat-transfer 

f h 

coefficient for the 1 sublump would be 


#4 = 


0.023 


r 0.4 

M 


D 0 ' 

1 


2 A i °- 


V - 8 


w °- 8 

l 


Hence the gross -lump heat flux becomes 


^t = C ht F W° l8 [T w -T b ] 


where the effective-geometry coefficient is 


C 


ht 


0.023 


l 



A ht,i 
D?’ 2 a?- 8 


and the fluid-properties coefficient is 


F = 


k 0.6 0.4 
°P 


M 


0.4 


* T b T b, av 


Effective-Geometry Coefficients 

The effective -geo me try coefficients are a, b, and c for the momentum equation 
and C ht for the energy equation. Each coefficient contains at least one summation, 
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which is performed over the length of fluid-lump. Consider, for instance, the volume 
summation 


v =i>i ax i 

i=l 

The l sublumps are chosen so that they correspond to lengths of component hardware 
with a constant cross-sectional area. The pump -discharge feedline, downstream from 
the main flow-control valve, is an example of a sublump. The cross-sectional area is 
uniform for an approximately 18-foot (5. 49-m) length. 

A curve of the total, accumulative volume for the axial -pump engine system from the 
pump exit to the core exit may be made (see fig. 28). The pump discharge is selected as 
the origin (the three -fluid-lump model begins at the pump discharge). The first linear 
segment of the curve represents the constant-area feedline from the pump exit to the 



0 3 6 9 12 15 18 

Distance, m 

Figure 28. - Accumulative volume of axial -pump system from pump 
exit to core exit 
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main flow -control valve. The remaining linear segments correspond to the other system 
components . 

The total volume in the liquid fluid-lump may be computed from 

l 

v tot,!< x s; )= y v i - 

i=f 

This represents the point on the volume curve in the figure at x = X g £. Clearly, the 
volumes of all three fluid-lumps are readily calculated from three points (X g ^ , X g g, and 
X ) on the accumulative volume curve; that is, 


X s I 



k s l 


v to M <v = 2 


x=0 


X, 


sg 


V tot,2p*8PV = E V x - V tot ,1 <*sZ> 
x=0 


L sg 


V tot,g^ X sg) V tot,^,2p,g (X ce) ^ V x 

x=0 


A similar argument may be made for each summation appearing in the effective - 
geometry coefficients. For each summation, accumulative curves may be derived, which 
represent the system from the pump discharge to the core exit. The value of each 
summation, therefore, is available at any point x (0<x<X ce ) . 

The effective -geometry coefficients are continuous functions of the positions of the 
two -phase fluid boundaries, X g ^ and X gg . In addition, the friction pressure -drop coef- 
ficient b is also a function of flow W. The coefficients may be written with the following 
functional notation. 
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Liquid fluid-lump: 


a l “ a ^ X sP 

b * = b < x 3- W i> 
C; = c(x a ) 


Two -phase fluid-lump: 


a 2p ** a ^ X sZ > V 






c 2p = c ^ X sZ’ V 


Gas fluid-lump: 


V b(X sg’V 
Cg = c(X sg ) 
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APPENDIX B 


TWO-PHASE EQUATION OF STATE 

The two -phase fluid is assumed to be a uniform, homogeneous mixture of liquid and 
gas phases. For such a mixture, an average specific enthalpy and an average specific 
volume are defined as 


h 2p X q h sg + t 1 h sZ 

(Bl) 

v 2 p = X , T sg. + [1 ’ VW 

(B2) 


where X is the mass quality, h g ^ = h g ^ (Pgp) saturated-liquid specific enthalpy, 

^sg = ^sg' P 2p^ is satu rated-gas specific enthalpy, v g ^ = v g ^ (Pgp) is the saturated- 
liquid specific volume, and v = v (P 9 ) is the saturated-gas specific volume. 
Eliminating mass quality X from equations (Bl) and (B2) results in 


v 2p = v sl 


y sg v s 1 

h + 

v sg v s 1 

^sg " h sz 

Sl + 

_ h sg - V 


2p 


(B3) 


Multiplying both sides of equation (B3) by Pg p yields 

P 2p v 2p = “< P 2P ) + ,3<P 2p )h 2p 


where 


and 


a = P, 


r 

_ 


|v - 

V sg ' v s l 

h sZ> 

ysz 

_ h sg - V 

= p„ 

v sg ■ v sl 


2p 

_ h sg " h sZ_ 



(B4) 


The variables a and /3 are well defined functions of the two -phase pressure Pg p and 
are shown in figure 29. 
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Equation (B4) represents the two -phase equation of state. Since 


P 


1 



equation (B4) may be rewritten in terms of the two -phase density, such that 

P 2p = p 2p[°' (P 2p> 

Given values for p 2p and h 2p , an algebraic balance of equation (B5) results in the 
proper value of P 2p . 
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APPENDIX C 


FRICTION FACTOR FOR TWO-PHASE PRESSURE-DROP EQUATION 

As discussed in the INTRODUCTION, the two-phase fluid was situated in the same 
general region of each engine system during the time periods of interest, namely, be- 
tween the nozzle-inlet manifold and the reflector exit. Because of the high metal tem- 
peratures of the reflector and the nozzle -coolant-tubes shell, the two-phase fluid-lump 
is treated as a high wall-to-bulk temperature ratio region. 

It is shown in reference 6 that the effect of high wall-to-bulk temperature ratios on 
the friction factor may be represented by the following correlation: 


f = 


Pi 


av 


where f| is determined by a modified Reynolds number in which the density and vis- 
cosity are evaluated at the film temperature (midway between the fluid and wall tem- 
peratures). Hence, for turbulent flow with Re < 200 000, an approximate expression 
for the two -phase friction factor is 
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APPENDIX D 


MATHEMATICAL REPRESENTATION OF ENGINE SYSTEMS 

Most of the variables are only functions of time. To simplify the equations, these 
variables are written without the function of time notation; that is, a typical variable 
y(t) is written simply y. However, some variables are either functions of both dis- 
tance and time, or are functions only of distance. These variables are typically written 
as y(x,t) and y(x), respectively. For clarity, this notation is used only in the final 
step of lengthy equation derivations. 


Tank 

The run tank is treated as a constant-pressure, liquid- hydro gen supply. Hence, 

Pj. = Constant 


Pump-Inlet Line 


The all-liquid pump-inlet line is a stationary, constant-geometry lump. The equa- 
tions derived in appendix A apply to this fluid-lump with stationary boundaries. Neglect- 
ing the momentum pressure drop yields 


3W • 


H = a „ 


at 


pi 


p . AZ . 
P x - P . - H 


pi 


n 


” b pi 


* ? 

W . 
Pi 


Pi 


For small variations in fluid enthalpy, 



the density p pi 


and the geometry coefficients 


a p ., b pi , and V pi are constants. 
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Pump 


The pump operation is assumed to be defined by a normalized, steady-state map. 
The maps utilized were derived from experimental data and are shown in figure 7. 
The set of defining equations for the pump is 


W 

P 


W . + W 7 
pi l 
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N = N(t) 


AP AP /W 
P PI P 
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P_ = P . + N 
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|AP_ 


N 


Pump Discharge to Core Exit 

The t hree-fluid-lump model . - The liquid fluid-lump is characterized by a station- 
ary upstream boundary at the pump discharge and by a variable downstream boundary 
at the saturated-liquid interface. The set of defining equations (appendix A) for the 
liquid fluid-lump is 


= Constant 


at 


= a ; [P. 


7 L po 
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‘ P 2pl ' b l — 




= b ( X sP^) 


*1 
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No mass -continuity equation is required since no mass-flow storage is assumed in this 
lump. An approximate value of nozzle-inlet-manifold pressure (for the high-speed 
analog) is 




where 



= f (X ni’ W P 

a 


The two -phase fluid-lump is bounded by two variable boundaries: the saturated- 
liquid boundary on the upstream side and the saturated-gas boundary on the downstream 
side. The boundary velocities Ug and Ug g g are approximated by analog deriva- 
tive circuits. The defining equations are 


~ [ p 2p, av V 2p] " W l “ p sl A sl U B, si " W 2p + p sg A sg U B, sg 
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b 2p “ b ( X sZ » X sg’ W 2p’ 


P f,2p 

P b, 2p y 


c 2p “ c ( X s£ » X sg) 


b av, 2p 


= 2p sl^sg_ 
p s l + p sg 


A two -phase fluid state relation is derived in appendix B. The resulting equation is 

P 2p = p 2p, av + ^ h 2p, avl 


where 


« = a ( p 2p) 


and 


p - «> V 

The gas fluid-lump is bounded by a variable upstream boundary at the saturated- 
gas interface and by a stationary downstream boundary at the core exit. The defining 
equations are 


JL[Y> VI = W„ - p A U_ - ,, 
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V b < x sg’V 

c g ° c(X sg> 

_ ^sgh nc 

av > g p + p 

H sg H nc 
P g = Rp g, av T b, g 

The system-heat-transfer model . - The operation of this part of the overall model 
is based on the conservation of energy in the components material. The system is con- 
sidered composed of many lumps. Each lump represents a fixed amount of mass. A 
typical mass -lump is represented in figure 30. 



Each lump is characterized by an average temperature T w j. Conduction from the 
outside surface to the inside surface is assumed negligible. Consider heat transfer 
between the metal and fluid only, and assume that the rate of heat loss by the material 
is balanced by the input heat flux to the fluid. 

From a heat balance for such a lump it follows that 


ST. 


Qi 


n t = - 


st 


M.c . 
1 p, 1 
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Let X. . be the distance from x = 0 (pump discharge) to one edge of the mass -lump 

1 , 1 

and X. 9 be the distance from x = 0 to the outer edge. Then 
b * 



A 2 Q' (x. 


t)dx 


This calculation will be performed for each lump in the spatial -distribution model. The 
inputs to this part of the model, therefore, are the individual mass -lump heat fluxes. 

The outputs are the average wall temperatures. 

The system -heat-transfer model is also programmed on the slow- speed analog. This 
part of the mathematical model considers the system from the nozzle-inlet manifold 
x = X n ^ to the thrust-nozzle chamber x = X n(; . A diagram of the mass lumping is shown 
as unshaded areas in figure 31. The convergent-divergent nozzle is treated as three 


x “ *ce 



mass-lumps. One mass-lump is assumed from the reflector inlet to the core exit. 

In the diagram the mass-lumps of the convergent-divergent nozzle are labeled 1, 2, 

and 3. Mass -lump 1 represents the outside nozzle -coolant-tubes material extending 

from the nozzle-inlet manifold x = X . to the nozzle-coolant-tubes throat x = X , . 

ni nt 

From the experimental data, the temperature of this tube material is about equal to the 
coolant fluid temperature. Hence, no heat flux is assumed in mass -lump 1. 

Mass -lump 2 represents the inside nozzle -coolant-tubes material extending the entire 
length of the tubes. This tube material is in contact with the coolant fluid on one side and 
with the thrust-nozzle hot gas on the other side. Based on experimental data, no heat-flux 
storage will be assumed. Hence, the heat flux gained from the hot gas is transferred 
entirely to the coolant fluid. The average metal temperature is assumed a constant: 
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T w, hg = Constant 


Mass-lump 3 represents the remaining outside nozzle -coolant-tubes material plus 
an external metal shell, or skirt. This mass-lump extends from the nozzle -coolant- 
tubes throat x = to the nozzle-coolant-tubes exit x = X ne . The defining equation is 


3T, 


at 



^ns c p, ns 


Geometry and mass data for the mass-lumps of the convergent -divergent nozzle are 
available in reference 7. 

The remainder of the system mass is grouped into a gross mass -lump 4. The initial 
system component composing mass -lump 4 is the reflector. The reflector is a source of 
high input heat flux. Therefore, it is anticipated that two -phase fluid will not penetrate 
much into the reflector. 

Most of the fluid in mass -lump 4 is assumed to be gas. The total heat flux lost by 
the mass in this lump is considered transferred to the gaseous fluid. Hence, the de- 
fining equation of mass -lump 4 is 


3T. 


w, n-ce 

at 


Q 

g 

^ri-ce c p, sg-ce 


The equations and operation of the slow-speed analog have been described. In 
review, a slow-speed solution is continually generated. Every second of computer 
operation produces 0.01 second of the solution. The pump speed and nozzle-inlet- 
manifold fluid enthalpy are programmed from experimental data. 

The slow-speed analog requires several inputs from the high-speed anaiog. These 
inputs are (1) positions of both two-phase saturated-fluid boundaries, (2) average fluid 
enthalpy for proper mass-flow storage in the two-phase fluid-lump, (3) average fluid 
temperature in the gas fluid-lump, (4) thrust-nozzle-chamber fluid temperature, (5) input 
heat flux from the nozzle-coolant-tubes shell, and (6) input heat flux to the gaseous fluid. 
A list of these inputs is X gl , X , h 2p , T b g , T fo nc , Q nS - and Q g . 

The slow-speed analog supplies several outputs to the high-speed analog. These 
outputs are (1) time, (2) nozzle-inlet-manifold pressure, (3) individual fluid-lump flows, 
(4) saturated-fluid densities, (5) average nozzle-coolant-tubes metal temperature for hot- 
gas heat transfer, (6) average metal temperature of the nozzle-coo«mt-tubes shell, and 
(7) average metal temperature of the region from the reflector inlet to the core exit. A 
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list of these outputs is t, P ni , W 2p , W g , p gZ , p gg , T w hg , ^ w>n8 , and T w> ri _ ce - 

The spatial-distribution model. - At some time t this part ot the model receives 
inputs from the slow-speed analog. Spatial-distributions of fluid enthalpy, pressure, 
density, and input heat flux are rapidly determined. From the distributions, certain 
quantities are calculated and sent to the slow-speed analog. 

At some later time (t + At), inputs are again received. Spatial distributions are 
again rapidly determined. And new outputs are calculated for use in the slow-speed 
analog. 

These new outputs represent updated values of the same quantities which were calcu- 
lated at time t. In the limit as At approached zero, the outputs are updated at each 
instant of time. 

This part of the model furnishes the three -fluid-lump model with the positions of the 
two -phase fluid boundaries in the system. It also supplies an average fluid enthalpy for 
the two -phase fluid-lump and two fluid temperatures for use in the gas fluid- lump. In 
return, the individual fluid-lump flows are received as inputs. 

Wall material temperatures are also received as inputs. The system -heat -transfer 
model furnishes these inputs, and receives, in return, evaluations of input heat flux for 
the individual mass -lumps. 

A complete cycle of operation will consist of four steps. The first step is the estab- 
lishment of all necessary initial conditions. These initial conditions are based on the 
inputs from the slow-speed analog. 

The second step is a continuous high-speed integration with respect to distance. 

The range of integration is from the nozzle-inlet manifold to the saturated-gas boundary. 
As a result of the integration, the following slow-speed-analog inputs are determined: 

(1) the two -phase saturated boundaries X g ^ and X gg , (2) the average two -phase fluid 
enthalpy h~ and (3) the input heat flux from shell of the nozzle-coolant tubes Q . 

“P j d v ns 

The third step of the operation cycle begins with the termination of the high-speed 

integration. The core-exit fluid temperature T fe ce is calculated from a heat balance 
on the gaseous fluid. The average values of the gas temperature T fe and input heat 
flux Q are also determined. At the end of this step, with all the spatial -distribution- 
model calculations completed, the new values are tracked, or sampled. The analog 
circuits (ref. 8) remain in the tracking mode only long enough to completely establish 
the new values. A typical output from a track-and-store circuit could take the form 
shown in figure 32. 

In the final step, information is transferred from the high-speed analog to the slow- 

speed analog, from the track-and-store circuits. The information remains stored until 

new values are calculated. The slow-speed analog inputs are transferred as soon as 

they are calculated. For example, one of the spatial distributions generated is that of 

fluid enthalpy h(x). During the high-speed integration, h continually increases from 

its initial value of h . to its final value of h . At some point (x = X.) h will pass 

ni sg st 
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Time 


Figure 32. - A typical output from a track-and -store circuit (greatly 
exaggerated). 


through the value of h = h g ^ . The instant that this occurs, a track-and-store circuit is 
triggered to track the value of x. The new position of the saturated-liquid fluid boundary 
x = X g £ , therefore, will be stored. 

The spatial- distribution model considers the system from the nozzle-inlet manifold 

(x = X .) to the core exit (x = X ). A diagram of the flow passages comprising the lump 
111 c e 

used in the simplified analog program is shown as the shaded areas in figure 31. 

The heat-transfer equations pertaining to the fluid in the nozzle -coolant tubes are 
based on the material presented in reference 7. Some of the expressions, of course, 
reflect assumptions already made for the slow-speed analog equations. 

The differential equations for the spatial-distribution model are developed in appen- 
dix C. The energy equation is 


Sh 

0X 



I 




I 


where 

<p(x, t) = <p t (h, P ni ) for h(x, t) < h sZ (x, t) 

= cp 2 p ( h , P 2p ) for h sl (x, t) < h(x, t) < h gg (x, t) 

subject to the initial condition 

h < X ni> « “ h ni<« * h ni 

The time derivative is approximated with an analog derivative circuit. 

The approximate pressure equation is 


ap 

3x 


(x, t) = - 


wDt 


w 


nA 


(x, t) 


The initial condition for this expression is supplied by the slow -speed analog and is 

P < x ni, t> - p rt<‘> = P „i 

The flow is assumed uniform in each fluid-lump of the slow-speed analog. Hence, 
the following discontinuous flow distribution is assumed: 


W(x, t) = W z 

for 

h(x, t) < h gZ (x, t) 

II 

to 

for 

h sZ (x,t) <h( Xj t) 

= w 

g 

for 

h(x, t) > h gg (x, t) 


The saturated-liquid fluid enthalpy h gZ is strictly a function of pressure* Because 
of a spatial pressure distribution, however, an equation may be written: 

h sl (P) = h c7 (X, t) 


si 


As in the slow-speed analog, the saturated-gas fluid enthalpy h is considered a con- 
stant. 


I 4 
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An expression for the effective wall shear r w for turbulent flow with Reynolds 
numbers less than 200 000 is (ref. 5) 


T w 


i fb P b u2 = °.02S 
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px f A 

LwdJ 




P b A' 


Bend and entrance effects may be included in an additional loss term . The effective wall 
shear becomes 
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For the case when variations in fluid viscosity are small. 


— (x,t) = - 

ax 


K fr (x)W 1 ‘ 8 (x,t) 
p(x, t) 


where 


K fr (x) 


0.023 ttDmJ- 2 
ngA^L 



In the liquid phase, the fluid density p is approximated by its saturated-liquid value. 
In the two -phase fluid, the state relation derived in appendix B is used. For the spatial - 
distribution model, therefore, the fluid density is 

P(x, t) =p sl (x, t) for h(x, t) < h gZ (x, t) 


P(x, t) 

a(P) + £(P)h(x, t) 


for h gZ (x, t) < h(x, t) < h gg (x, t) 
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The input heat flux per unit length Q'(x) takes on three basic forms. The describ- 
ing equation is 


Q’(x, t) = Q^g(x, t) for X ni < x < X nt 

= Qhg (x > + ^is (x ’ *> for X nt " x - X ri 
= Qp(x,t) for X ri <x<X re 

where Qjjg( x ) t) is the input heat flux per unit length from the hot gas on the thrust side 

of nozzle, Q* (x,t) is the input heat flux per unit length from the nozzle -coolant -tubes 
ns 

shell, and Q^(x,t) is the input heat flux per unit length from the reflector passages. 
The equation for hot-gas input heat flux is 

^hg “ ^hg^ht, hg t^b, nc “ ^w, hg-1 

where Aj^ is the hot-gas heat-transfer area per unit length. The hot-gas heat- 
transfer coefficient will be similar to that used in reference 7. The expression 

for the coefficient is 


^hg" C hg~ Re f Pr f 


where 


Re^ = 


WD 
Mf A 


The input heat-flux equation, therefore, may be written as 



where 




nO. 8 
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ht, hg 


St^O. 4-, 


F hg = [kM Pr ] T ^_ T 


b f, hg 


and 


T + T 

t _ w, hg b, nc 
f ,hg 2 

The fluid-properties coefficient F^g is evaluated at an average, constant value of 
pressure. The geometry coefficient C^gCx) was obtained from reference 7. 

The nozzle-coolant -tubes wall between the hot-gas and the coolant is assumed to be at 
a constant temperature T w ^ . By assuming no heat-flux storage in the wall Q^g may 
be calculated and applied to the coolant side of the nozzle. Hence, no complicated two- 
phase boiling heat-transfer coefficient is required. 

The input heat flux from the nozzle-coolant -tubes shell Q ng however, requires the 
use of a boiling heat -transfer coefficient. The convective film boiling heat-transfer 
coefficient from reference 9 is used. The correlation is for liquid hydrogen in the pres- 
sure range 25 to 175 psia (172 to 2107 kN/m ). In terms of the Martinelli parameter 
Xy., the correlation is 


NUexp > f = + 0. 15 

Nu calc, fm °- 7 + 2 ’4 
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p fm p f p s l 


T f = 


T w + T b 
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x tt " 
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The input heat flux per unit length from the nozzle-coolant-tubes shell is 

^ns “ ^ns^ht, ns !- T w, ns " T b, ns^ 

where ng is the heat-transfer area per unit length. The bulk fluid temperature 
given by 

^b,ns = T sat^ 

Combining the last six equations results in the following expressions 
^s^*) = K ns< x > F ns' x > t )' S ' 0 ’ 8 ( x " l )t T w,ns ‘ T b, rj <*> « 
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and 


T + T 

ry, „ w, ns b, ns 

f no 

2 


f, ns 


The fluid-properties coefficient F is evaluated at an average, constant value of 

Ufa 7 

film temperature T f __ at several values of pressure. The final value of F is 
obtained from linear interpolation between the constant -pressure values. The wall tem- 
perature T is an input from the slow-speed analog and is assumed uniform through- 

W y US D 

out the nozzle -coolant-tubes shell. 

The overall nozzle shell input heat flux is the sum of <^ g for the entire 

length. Hence 


Q =/ 

tis 


X ri 


X, 


nt 


^ns(x, t)dx 


When x = X pi during the high-speed integration, this value is stored for use in the 
slow-speed analog. 

The final input heat flux per unit length to be defined is Q^., the value for the two- 
phase fluid in the reflector passages. It is assumed that the input heat flux per unit 
length of the reflector -inlet manifold is approximately the same as that in the reflector 
passages. Two-phase film boiling is assumed. The defining equations are similar to 
those for the nozzle-coolant -tubes shell; that is, 


Q’ (x, t) = K (x)F (x, t)W°‘ 8 (x, t)[T 


sat 


](x, t) 


where 


K r = 


0.023 

D 


0 . 8 .» 
A ht,r 


F r = 


[kM-°* 8 Pr 0 * 4 ] fj 

r ^ 

p sZ p f 

0.8 

L 

1. 105 + 0. 36 x tt " 

0. 8 
p b 

Vs l + I 1 - X q ]P t 

L J 

r 

0.7 + 2. 4 Xfl. 
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T _i_ t 

rp _ w, r sat 
f,r „ 


and 


T = T - AT 

w, r w, ri-ce r, 2p 

The fluid properties coefficient F r is evaluated at an average, constant film tem- 
perature T f at several values of pressure. The final value of F is obtained from 
a linear interpolation between the constant-pressure values. 

As discussed earlier, two-phase fluid is found just inside the reflector. The tem- 
perature of the reflector material in contact with the two -phase fluid is T . The 

w, r 

T w r is assumed to be lower than the average wall temperature T w r ^_ ce by a fixed 
amount AT r gp- The T w r i_ ce i s the temperature of mass-lump 4 in the slow-speed 
analog. This is a good approximation as shown in figure 16. The difference between the 
predicted T w r j_ ce and the experimental reflector -inlet fluid temperature tends to 
remain constant. 

The high-speed integration ends at the saturated-gas boundary x = X . The re- 

sg 

mainder of the system contains gaseous fluid. A general equation for the overall input 
heat flux to the fluid is 


^g ~ ^g A ht, gI T w, ri-ce T b, gl 

The conventional correlation of reference 5 for gaseous hydrogen flow is used; that is, 

Nu b = 0.023 Re b - 8 Pr b * 4 

Combining the two previous equations and then grouping terms yields the following 
relation for the overall, gaseous -fluid, input heat flux: 

= C ht, g (X sg >F g'' , g 8 [ T w, ri -ee " T b, gl 
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where 


ce 




0.023 A 


ht 


x=X 


D °. 2 a 0 . 8 


sg 


F g = [> c (‘' 0 ' 8 Pr 0 - 4 ] T T 
5 b b, g, av 


b, g, av 


_ T sat + T b, ce 


and 


T sat “ T sat^ P g^ 


The derivation of the effective -geometry coefficient is presented in appen- 

dix A. The fluid-properties coefficient is assumed a constant; its value is deter- 
mined from average, constant experimental values of T^ ay and gas pressure. 

The input heat flux per unit length Q* is required. An approximate expression is 


Q 

Q’ -—3 1 

g X 


V 

ce sg 


The only remaining term to be evaluated is the thrust-nozzle -chamber fluid tem- 
perature T^ . The core-exit fluid temperature may be obtained from an energy bal- 
ance on the gaseous fluid. A simplified finite-difference form of the general energy 
equation found in appendix E is used. Neglecting the time derivatives and assuming 
uniform flow, 


'■M, 


Ql 
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This equation may be rewritten as 


W c 
S P, g 


T - T 
b, ce sat 

- X_ 
ce sg _ 



The thrust-nozzle-chamber fluid temperature may be approximated by the value of the 
core -exit fluid temperature; that is, 


b, nc 


T k = T o , + 
b, ce sat 


*g c p,g 


The specific heat Cp is evaluated at an average, constant value of fluid temperature. 

The description of the high-speed analog is nearly complete. The necessary 
equations for the high-speed solution have been presented. The averaging procedure for 
h 9r _ o,r> the average two-phase enthalpy, however, remains to be described. Assume 
that this averaging is not required in the gas fluid-lump . 

To obtain proper mass-flow storage in the gross two -phase fluid-lump an average 
value of the thermodynamic property [Bp/BP]^ is required. This phenomenon was dis- 
cussed earlier in appendix A. Constant -pressure curves of [0p/0P] ^ against fluid 
enthalpy are derivable from hydrogen properties tables (ref. 3). 

Spatial distributions of fluid enthalpy h and fluid pressure P are continuously 
generated every cycle. By linearily interpolating on pressure, a continuous distribution 
of [9p/3P] k is derived from the enthalpy distribution. Hence, 



is known. 

The average value of this thermodynamic property for the two-phase fluid-lump 
would be 




Therefore, the average two -phase fluid enthalpy for proper mass-flow storage is 



This required property of enthalpy as a function of pressure and [dp/d P] b is derivable 
from the curves in reference 3. 

The high-speed analog requires several inputs from the slow-speed analog. The 
inputs are (1) time, (2) nozzle-inlet-manifold pressure, (3) individual fluid-lump flows, 
(4) saturated-fluid densities, (5) average nozzle-coolant -tubes metal temperature for 
hot-gas heat transfer, (6) average metal temperature of the nozzle-coolant-tubes shell, 
and (7) average metal temperature of the reflector -to -core-exit region. A list of these 

inputs is t, P n j, W z , W 2 p, W g , P gZ > b S g> ^w, hg’ ^w.ns’ anc * ^w.ri-ce' 

The high-speed analog supplies several outputs to the slow-speed analog. These 
outputs are (1) the positions of both two-phase saturated fluid boundaries, (2) the average 
fluid enthalpy for proper mass-flow storage in the two-phase fluid-lump, (3) the average 
gas -fluid -lump temperature, (4) the thrust-nozzle-chamber fluid temperature, (5) the 
input heat flux from the nozzle-coolant-tubes shell, and (6) the input heat flux to the 
gaseous fluid. A list of these outputs is X gZ , X gg , h 2p ay , T^ g , T b nc , Q ns , and 

Qg- 

Each cycle of operation begins with the establishment of initial conditions. These 
conditions are based on inputs from the slow-speed analog. During the subsequent high- 
speed x-integration, updated values for X gZ , X gg , h 2p Q ns are formed and stored 
in analog track-and-store circuits. And finally, new values for T b , T b , and Q g 
are calculated and sent to the slow-speed analog. 


Thrust Nozzle 

The analog model is terminated by a thrust -nozzle lump. The thrust nozzle is 
assumed to be choked. The defining equations are 


^nc _ W g - W nc 
9t V nc 


w = 

nc 


K „c P nc 




nc 
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and 


P„c = Rp nc T b,nc 
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APPENDIX E 


SPATIAL-DISTRIBUTION MODEL EQUATIONS 

A set of general, one -dimensional differential equations based on the conservation of 
mass, momentum, and energy, respectively, in circular passages is as follows: 

+ 1 9W = o 

at a 9x 

ffDgr w 
A 

— [ph-j’P] +IJL[Whl 
at A 0x A 

Combining the mass and energy equation yields 

p a h_j’ 9£ + W0h = Q^ 
at at A 0x A 

Solving for the spatial gradient yields 


1 0W 0 


A at ax 


ngP + 




pA‘ 


0h 

0X 


Q’ + A 


T* 0P _ 0h 
J p — 

at at 
w 


Neglecting mass-flow storage results in 


0h 

at 


PS 


^h 0P 
3PJp 3* 


since 


0p _ _I 9W w Q 

at a ax 
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The enthalpy spatial gradient for fluid-lump, therefore, becomes 


where 


Q’(x, t) + A (x)<o(x, fc)— (x, t) 

— (x,t) = — d ± 

3x W(x, t) 


<p(x, t ) = j' - p(x, t) py (x, t> 

L3pJ p 


and, to a first approximation, 


— (x, t) 


aT3 3P~ 

IE (t) = 

3t 3t 


For the pressure ranges in this study the pressure dependence of the saturated-gas 
fluid enthalpy may be neglected (ref. 3) . The pressure dependence of the saturated- 
liquid enthalpy, however, is important. The pressure-drop equation is integrated with 
respect to x in repetitive operation solely to supply this pressure dependence. There- 
fore, to generate an approximate two-phase pressure distribution and, at the same time, 
minimize the necessary amount of analog computing equipment, only friction pressure 
drop is considered; that is, 


^<x,t) = - 
3x 


ttDt 


w 


nA 


(x,t) 
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APPENDIX F 


SYMBOLS 

A flow area, ft 2 ; m 2 

t 2 2 

A ht heat -transfer area per unit length, ft /ft; m /m 

a effective-geometry coefficient, (lbm)(in. 2)/(ibf)(sec 2 ); (kg)(m 2 )/(N)(sec 2 ) 

-3 -3 

b effective -geometry coefficient, ft ; m 

C geometry coefficient 

-3 -3 

c effective-geometry coefficient, ft ; m 

c specific heat of hydrogen at constant pressure, Btu/(lbm)(°R); J/(g)(K) 

P 

D diameter, ft; m 

F fluid properties coefficient, Btu/(lbm°' ®)(ft®* 2 )(sec 0, 2 )(°R); kJ/(kg°* 8 )(m°* 2 ) 

(sec 0, 2 )(K) 

f coefficient of friction 

g gravitational acceleration, 32. 174 (lbm)(ft)/(lbf)(sec 2 ); 1 (kg)(m)/(N)(sec 2 ) 

heat-transfer coefficient 
h specific enthalpy, Btu/lbm; J/g 

J f conversion factor, (144 in. 2 /ft 2 )/(778 ft-lbf/Btu) 

K heat-transfer geometry coefficient 

Kfogg pressure loss coefficient for bends and entrance effects 
k thermal conductivity of hydrogen, Btu/(ft)(sec)(°R); kJ/(m)(sec)(K) 

L length, ft; m 

l number of sublumps 

M mass, lbm; kg 

m pump -map slope, (lbf)(sec)/(lbm)(in. )(rpm); (N)(sec)/(kg)(m )(rpm) 

N speed, rpm 

Nu Nusselt number 

N u ca ic calculated Nusselt number for two-phase flow 
Nu experimental Nusselt number for two-phase flow 

n conversion factor, 144 in. /ft 
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II ‘ ' 

n 

P pressure, psia; N/m 

Pr Prandtl number 

Q input heat flux, Btu/sec; kJ/sec 

Q' input heat flux per unit length, Btu/(ft)(sec); kJ/(m)(sec) 

R specific gas law constant for hydrogen, 5.326 ftV(in. ^)(°R); 420.8 m/K 
Re Reynolds number 

T temperature, °R; K 

t time, sec 

U velocity, ft/sec; m/sec 

V volume, ft**; 

v specific volume, ft /lbm; m /kg 

W mass -flow rate, lbm/sec; kg/sec 

X distance (finite length), ft; m 

X mass quality 

T. 

x distance (independent variable), ft; m 

Z height, ft; m 

a two-phase state -equation coefficient, (lbf)(ft )/(lbm)(in. ) abs; (N)(m)/kg 

/3 two-phase state -equation coefficient, (lbf)(ft )/(Btu)(in. ) abs; (N)(m)/kg 

p viscosity, (lbm)/(ft)(sec); kg/ (m) (sec) 

p fluid density, lbm/ft^; kg/m^ 

P av average (linear) fluid density, lbm/ft ; kg/m 
r effective shear, lbf/ft^; (N/m^) 

cp heat-transfer function, (Btu) (in . ^) / (lbf) (ft^) 

X^. Martinelli two -phase parameter 
Subscripts : 
av average 

B boundary 

b bulk 

ce core exit 
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ci 

core inlet 

e 

error 

eff 

effective 

f 

film 

fm 

film mean 

fr 

friction 

g 

gas 

gr 

gravity 

h 

enthalpy 

hg 

hot gas 

ht 

heat transfer 

i 

sublump number 

in 

inlet of flow stage 

iner 

inertia 

l 

liquid 

loss 

loss 

mfcv 

main flow-control valve 

mom 

momentum 

n 

nozzle 

nc 

thrust -nozzle chamber 

ne 

nozzle-coolant-tubes exit 

ni 

nozzle -coolant-tubes inlet 

ns 

nozzle -coolant-tubes shell 

nt 

nozzle-coolant-tubes throat 

out 

outlet of flow stage 

P 

pressure 

P 

pump 

Pi 

pump inlet 

PO 

pump outlet 

r 

reflector 
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re reflector exit 

ri reflector inlet 

sat saturated 

se spider exit 
sg saturated gas 

sZ saturated liquid 

t tank 

tot total 

w wall 

wb wall to bulk 
x distance 

2p two -phase 
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